An entropy perspective on the thermal crossover in a fermionic Hubbard chain 



m 

o 

(N 

Oh 
< 

0^ 



I 



I 

O 

o 



> 

o 

en 
in 

^' 

o 
en 

> 



Lars Bonnes,'* Hannes Pichler,'-^ and Andreas M. Lauchli' 

^ Institute for Theoretical Physics, University of Innsbruck, A-6020 Innsbruck, Austria. 
Institute for Quantum Optics and Quantum Information of the Austrian Academy of Sciences, A-6020 Innsbruck, Austria. 

(Dated: April 22, 2013) 

We study the Renyi entropy in the finite temperature crossover regime of a Hubbard chain using quantum 
Monte Carlo. The ground state entropy has characteristic features such as a logarithmic divergence with block 
size and 2fcF oscillations that are a hallmark of its Luttinger liquid nature. The interplay between the (extensive) 
thermal entropy and the ground state features is studied and we analyze the temperature induced decay of the 
amplitude of the oscillations as well as the scaling of the purity. Furthermore, we show how the spin and charge 
velocities can be extracted from the temperature dependence of the Renyi entropy, bridging our findings to recent 
experimental proposals on how to implement the measurement of Renyi entropies in cold atom system. Studying 
the Renyi mutual information, we also demonstrate how constraints such as particle number conservation can 
induce persistent correlations visible in the mutual information even at high temperature. 



I. INTRODUCTION 



Entanglement in quantum systems is a manifestation of in- 
herent quantum correlations between different degrees of free- 
dom'. Over the last decades, quantum information theory 
lead to a deeper understanding of how to characterize and 
quantify entanglement in few and many body systems with 
a wide range of applications in condensed matter and cold 
atomic systems. One of the most striking features for ground 
states for a wide class of physics Hamiltonians is the area law^ 
for the entanglement entropy stating that the entanglement be- 
tween two subsystems is determined by their boundary - alike 
the entropy of a black hole being proportional to the size of its 
horizon-' - and not the their volume. This has striking conse- 
quences for gapped one dimensional quantum systems where 
the entanglement is independent of block size"* consequently 
resulting in the success of the density matrix renormalization 
group^"^ (DMRG) and its descendants. Measuring entangle- 
ment entropies in real-life experiments, however, is quite chal- 
lenging since the Renyi entropies, for instance, are non-linear 
functionals of the reduced density matrix. The latter has actu- 
ally been measured using state tomography in trapped ions^-'" 
and superconducting qubits" for few particles but scaling this 
approach to systems with possibly hundreds or thousands of 
degrees of freedoms is challenging'^''''. 

Cold atom systems are promising candidates to implement 
scalable entanglement measurements as they provide the pos- 
sibility to engineer interacting quantum (lattice) systems and, 
as recent theoretical proposal show'** ''^, directly access the 
Renyi entropies of the system. In particular, as put forward 
by Daley et al. '^, the nth Renyi entropy for a bosonic lattice 
system can be infeiTed from the expectation value of the swap 
operator by performing correlated site-resolved density mea- 
surements on n different realizations of the systems. These 
measurements can be implemented by using a quantum gas 
microscope, as it has already been realized by different ex- 
perimental groups'^'^", that allows for a single-site resolved 
detection of local occupations. In a recent paper"*, we gener- 
alized this scheme to the case of fermionic atoms where spe- 
cial care has to be taken of phase factors from the braiding 
of two identical fermions. These new experimental tools can 



not only test theoretical concept such as dynamical entropy 
generation in quantum quenches^' or the conformal structure 
of critical one-dimensional systems but also provide insight 
in outstanding theoretical questing as corrections to the area 
law in (2 + 1) dimensions^^"^^ or topologically ordered sys- 
tems-^"-^'. 

Whereas conformal field theory is a powerful tool to obtain 
the entanglement properties of many-body systems in one di- 
mension, the crossover from a highly entangled many-body 
wave function to a thermal state as temperature is increased is 
much less explored. With the prospect of future experiments 
that will operate in this intermediate regime a deeper under- 
standing of the interplay between entropy and thermal entropy 
is of great interest. 

To address these questions, we use recent finite-temperature 
quantum Monte Carlo methods complemented by ground- 
state DMRG to access the Renyi entropies in a fermionic Hub- 
bard chain. We can track the ground-state features such as 
the logarithmic divergence of the entanglement entropy and 
the spatial 2k-p oscillations over a wide range of temperatures 
and demonstrate how the low energy Luttinger liquid can be 
probed beyond the determination of the central charge by ex- 
tracting the spin and charge velocities from the temperature 
behavior Renyi entropies. We also study the mutual informa- 
tion and show how constraints such as particle number con- 
servation affect the fluctuations that lead to finite mutual in- 
formation of a bipartite splitting even at infinite temperatures. 
From the perspective of prospective experiments, these results 
provide valuable insight on the temperature regimes where ef- 
fective ground state entanglement will be visible and what in- 
formation about the true ground state can be inferred from 
thermal Renyi entropies. 

This paper is structured as follows. In Section 11 A and 
II B, the one-dimensional Hubbard model and previous re- 
sults for the ground-state entanglement entropy and thermo- 
dynamic properties of the underlying Luttinger liquid are in- 
troduced. An additional measure for mixed states, namely the 
Renyi mutual information, is introduced in Section II C fol- 
lowed by a short discussion of the QMC methods employed to 
extract the Renyi entropies for thermal states in Section II D. 
The simulation results are presented in Section III where we 
show the finite-temperature crossover in Section III A and the 



extraction of the Luttinger velocities in Section IIIB. Finally, 
Section III C discusses the mutual information for constrained 
(canonical) as well as unconstrained system. 



II. MODEL AND METHODS 

A. One dimensional Hubbard model 

We consider spinful fermions on a chain of length L with 
open boundary conditions. One of the simplest interacting 
models is the Hubbard model with nearest-neighbor hopping 
and on-site repulsive charge interaction, [/ > 0, reading 
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The fermionic operators c\ ^ and c, ^^ create or annihilate a 
particle with spin a on site i and rii^a — 0,1 counts the num- 
ber of spin-(T fermions on that lattice position. Furthermore, 
we will only consider the spin balanced sector, i.e. n^ = n^ 
(n^ = l/LY,ini^a). 

The Hubbard model is exactly solvable via the Bethe 
ansatz^^'^^"* and exhibits a metallic ground state for generic fill- 
ings n = n^ + ni^l that is described by a two-channel Lut- 
tinger liquid of the spin and charge degrees of freedom with 



velocities Va and v, 



33-35 



At half filling, on the other hand. 



a charge gap opens for C/ > whereas the spin degrees of 
freedom remain critical - the system forms a Mott insulator 
Although the charge degrees of freedom are frozen out the dy- 
namics of the spins give rise to an effective antiferromagnetic 
superexchange coupling J = At^ /U of adjacent spins. 



B. Entropy 

The nth Renyi entropy of a subsystem A containing I a lat- 
tice sites is defined in terms of the moments of the reduced 
density matrix pA — tr^p as 
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logtrp^. 
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Here A denotes the complement of A and p is a not neces- 
sarily pure density matrix. For ?i — > 1 one recovers the von 
Neumann entropy 5'vn = ^tr[p^logpyi]. 

Insight into the entanglement properties and the logarithmic 
corrections to the area law for ground states of the Hubbard 
chain is provided by conformal field theory.^*-^^ The lead- 
ing contribution for a bipartite splitting of a chain with open 
boundary conditions reads"*^^" 



S^'^{A) = - ( 1 



1 



log2Z^ + const., 



(3) 



where l'^ — L / it sni{'Kl a / L) is the chord distance of block 
A. The prefactor of the logarithmic divergence is sensitive to 
the central charge and readily distinguishes between the the 



c—1 Mott insulator^"* and the semi-direct product of the two 
Luttinger liquids governing the metallic phase resulting in an 
effective central charge of c = 2^^'*'''^. In the case of open 
boundary conditions considered here, non-leading lattice ef- 
fects will lead to 2k-p oscillations even for the von Neumann 
entropy - here fcp is the Fermi momentum - whose decay 
with Z^ is governed by the Luttinger parameters"***'^""^^. Fur- 
thermore, there exist other corrections to Equation (3) such as 
shell filling effects."*^ 

Conformal field theory can also access the entropy of sys- 
tems at finite temperature and the von Neumann entropy for 
the Hubbard chain in the thermodynamic limit is given as^*'*'* 
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where Vc and Vg denote charge and spin velocities. This for- 
mula interpolates between the logarithmic growth with block 
size for T — > and the volume law scaling S'vn oc I a of the 
thermal entropy. A general expression for the Renyi entropies 
for a system of finite size and finite temperature has, how- 
ever, only been derived for massless Dirac fermions**^. Also, 
Ref 46 discusses the finite-temperature crossover in the finite- 
size XX-chain using correlation matrix methods. 

To make a connection to the thermodynamics of the Hub- 
bard model, we reexpress the thermal Renyi entropy for a 
mixed state p = exp{—H/T)/Z of the entire system in terms 
of the Helmholtz free energy F ^ —T log Z, namely 



C(r) = 



F{T/n) - F{T) 
T(l - n) 



(5) 



Focusing on the Mott insulating phase where the charge de- 
grees of freedom are gapped out, the low-temperature spe- 
cific heat for a Ia site system is just Cy = 7rr?^/(3us) 
where Vg = Ji{nU/{8t))/Jo{TrU/{8t)) is the spin veloc- 
ity '*^^'' and Jn{x) denote the modified Bessel functions of 
first kind. From these considerations one infers that the ther- 
mal n — 2 Renyi entropy in the T — > limit takes the form 



^f(r) = ^TZ^/(4i;,). 



C. Mutual information 



(6) 



The Renyi entropies for thermal states do not obey an 
area law due to the extensive contribution of the thermody- 
namic entropy and in general provide no information about the 
amount of information one subsystem has about the other. Un- 
like entanglement entropy, thermal entropy is merely a mea- 
sure of our lack of knowledge about the micro state than of ac- 
tual correlations. The mutual information /„, defined in terms 
of Renyi entropies as 



/„(A, B) = S,M) + SJB) - S,,{A U B), 



(7) 



on the other hand is sensitive towards correlations between 
two blocks^" A and B, i.e. not only two-point correlation 



functions that are used characterize possible (quasi) long- 
range order but all information one subsystem has about the 
other. For the bipartition of a pure state one simply recov- 
ers the Renyi entropy, i.e. In{A, A) ~ 2Sn{A), that is con- 
strained by the area law. In addition to that, it has also been 
shown that the mutual information obeys an area law even for 
thermal states^'^'"^^. 

The mutual information for thermal (classical) states has 
attracted a lot of interest recently with respect to phase tran- 
sitions in classical^''''''' and quantum spin models^^. Further- 
more, the /„ {A, A) is closely related to the excess entropy 
used in classical complexity theory as a measure for spatial 
memory of an infinite chain''*'. Thus, being able to access I2 in 
the experiment also opens new directions of study apart from 
aspect of (quantum) entanglement such as the loss of mutual 
information in a quantum quench. ^^ 



D. Numerical Methods 

The Renyi entropy can be accessed in QMC by using the 
idea the the 71th moment of the reduced density matrix can 
be written in terms of the ratio of two partition functions 
as trp^ = Zn{A)/Z'\^^ Here, Z ^ trexp(-iy/T) is the 
usual partition function and Z„ (A) lives on the n-sheeted Rie- 
mann surface with cuts along subsystem A (see e.g. Ref. 21 
and 36 for an illustration). Several schemes have been de- 
vised and applied to measure the partition function ratio us- 
ing classical Monte Carlo^*^ ''", valence-bond projector Monte 
Carlo*"', stochastic series expansion algorithms^'^'^^-*"^ or path- 
integral Monte Carlo^**. Here, we follow the ideas of Ref. 28 
of using a global topology update Zn{A) o Z" that al- 
lows us to directly sample Z„(A)/Z" by simply counting 
how long in terms of Monte Carlo time the system is in 
one or the other world line topology. The QMC simula- 
tions of the Hubbard chain at fixed filling are performed using 
a sign-problem free directed loop algorithm in the stochas- 
tic series expansion (SSE) framework''^^"*'^ and projecting the 
measurements into the subspace of fixed density. To avoid 
sampling problems due to small matching probabilities of 
the world lines across the cuts of Z„, we use the incre- 
ment trick^^ that splits the calculation for a block of size 
I A into the calculation of smaller blocks. This can be done 
by expanding the partition function ratio as Z„(Z^)/Z" = 

[Zr,{lA)/ZrXlA-5l)]..\Z.a{25l)/Z,-,{6l)][Zr,{5l)/Zr,m2.nA 

identifying Z„(0) = Z" . Each quotient can be simulated in- 
dependently and (5/^5- similar to the results by Humeniuk 
and Roscilde^** - is found to be a good tradeoff between the 
number of increments and simulation efficiency. 

It is noteworthy that Sn{A U B) can be readily obtained for 
arbitrary intervals A and B in QMC since placing the cut(s) 
in Zn is arbitrary whereas the extraction of l2{A, B) for arbi- 
trary blocks in DMRG is, in general, more involved. 

In addition to our finite-T QMC simulations DMRG^ '^ is 
used to access the ground state properties of the Hubbard 
chain. The results shown are well converged using up to 3000 
basis states for the largest systems of i = 64 sites. 




FIG. 1. (Color online) n = 2 Renyi entropy profile for the Mott 
insulator at U /t = 1 for a chain of length L = 48. 
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FIG. 2. (Color online) Purity P as a function of TjJ {T/t) for the 
half-(quarter-) filled chain with L = 48 sites. The error bars are 
omitted for clarity of the figure. Inset: System size scaling of the 
purity for the Mott insulator with U /t = 8.5. 



m. FINITE-TEMPERATURE CROSSOVER REGIME 
A. Entropy profiles and the purity 

Although the Luttinger liquid ground state of the Hubbard 
model is gapless in thermodynamic limit, a finite spatial extent 
and thus an infra-red cut-off leads to a gap to the first excited 
state that closes like 1/L for large system sizes^^. Conse- 
quently, the ground state can effectively be probed even at fi- 
nite temperatures as long as the temperature is well below the 
finite-size gap. For larger temperatures, the low lying linearly 
dispersing excitations will be populated leading to a finite cor- 
relation length in the system as well as an extensive thermal 
contribution to the entropy. 

First, let us focus on the Mott insulating case at half fill- 
ing. Superexchange leads to an effective description in terms 
of a spin model and the ground state exhibits 2kp correlations 
in the spin channel that are also reflected in the behavior of 



the Renyi entropy "*'. Figure 1 exemplifies the n — 2 Renyi 
entropy profiles for U/t = 1 (see Figure 4 for U/t — 6). 
For T/t < 0.02 the data for L ^ AS, agrees perfectly with 
the ground state DMRG results and the system is effectively 
at zero temperature. A second feature in the low-temperature 
crossover regime is the reduction of the amplitude of the 2/cf 
parity oscillation in the Renyi entropy. We use a heuristic fit- 
ting ansatz where the thermal and quantum contribution to the 
entropy are simply added (including subleading oscillations). 
This fit shows that the amplitude of the oscillations is sup- 
pressed exponentially in T/t compatible with the presence of 
a finite thermal correlation length. Above T/t > 0.2, how- 
ever, the oscillations are absent within our numerical accuracy 
but at the very boundary of the chain and the entropy increase 
is linear with blocks size. The system is then in a regime dom- 
inated by thermal fluctuations although the temperatures for 
the data shown in Figure 1 are still below the superexchange 
scale. 

Upon increasing the temperature, the entropy profile be- 
comes asymmetric across the center of the system and the 
purity P — exp[— S'2(£)] starts to decrease. It can be di- 
rectly read of from the entropy profiles and we summarize its 
temperature dependence for different values of the interaction 
and i = 48 in Figure 2. For larger values of U/t we find 
that the purity for T/J < 0.5 is independent of the interac- 
tion but shows a common decrease that is exponential in T/J 
above the finite-size gap. This can be directly traced back 
to the effective spin physics in the low energy sector that de- 
pends on the spin velocity that is approximately proportional 
to J for large [/. At T « J, P begins to flatten and the sys- 
tems enters the regime where the spin channel is disordered 
but the charge channel remains frozen since the temperature 
is still small compared to the interaction energy. For U/t — 1, 
however, things look differently because i) the spin velocity 
is about a factor of three larger as compared to U/t = 9 and 
ii) the extent of the system is still too small compared to the 
charge correlation length such that the system in only on the 
verge of being a c = 1 Mott insulator. The system size de- 
pendence of the purity is exemplified in the inset of Figure 2 
where the extensivity of the thermal entropy directly leads to 
an exponential suppression of P as a function of L. A de- 
tailed analysis of the connection between 5*2 (P) and Vs will 
be presented in Section III B below. 

The 5*2 profile for the metallic case is given in Figure 3 
where we focus on the quarter-filled case. The picture is sim- 
ilar compared to the Mott insulator regime. For system sizes 
larger than 48, however, we do not reach the (effective) ground 
state even for T/t — 0.02 opposed to the U/t — 1 data shown 
in Figure 1 . 

If the two velocities separate and, in particular, the charge 
velocity is much larger compared to Vg, it is possible to ac- 
cess the spin-incoherent Luttinger liquid regime*'^"*''^ where 
the spin channel is thermally disordered but the charge de- 
grees of freedom can effectively be described by their low- 
energy field theory. These effects will be reflected in a change 
of slope of the purity once the spin channel is completely dis- 
ordered. 
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FIG. 3. (Color online) n = 2 Renyi entropy profile for the quarter 
filled Hubbard chain at quarter filling and U/t — 4. The upper and 
lower panel correspond to 48 and 64 lattice sites respectively. 



n 


U/t 


QMC, S{L/2) 


QMC, S{L) 


exact^^-*^^ 


1 


6 


0.92 ±0.05 


0.937 ±0.01 


0.9257 


1 


8.5 


0.702 ±0.005 


0.720 ±0.006 


0.6929 


1 


9 


0.63 ±0.08 


0.64 ±0.02 


0.6588 


1/2 


4 


0.59 ±0.01 


0.6 ±0.1 


0.5961 



TABLE I. Vs and l/[l/vs + l/«c] extracted from the Renyi entropies 
for the Mott insulator and the quarter-filled metallic case from QMC 
simulations using either the half-system or the full-system entropy 
compared to the exact results. 



B. Extracting Vs and Vc 

The low-energy thermodynamics of a Luttinger liquid is 
governed by the spin and charge velocity and the prefactor of 
the volume law for the Renyi entropies encodes this informa- 
tion, as seen in Equation (6). Since we do not know an exact 
expression for the Renyi profiles at finite temperature, we ex- 
tract the velocities by either considering the entropy of the full 
system and the half system. The latter quantity is most insen- 
sitive towards boundary effects and logarithmic contributions 
from the quantum entanglement. 

Figure 4 shows the Renyi entropies for n = 1, U/t = 6 
and T/t = 0.06 for different system sizes. One can see that 
52 for large block sizes eventually converges towards the ther- 
mal volume law as it is expected from Equation (4) and (6). 
Fixing the lattice position L/2, S2/T is independent of tem- 
perature over a wide range of system sizes, as shown in the 
inset of Figure 4. Using a linear fit for system sizes L > 32, 
we obtain the spin velocities from our ansatz in Equation (6) 
tabulated in Table I. A similar procedure is used to extract the 
velocities from the full system entropy (see inset in Figure 2). 
For the quarter filled system, we simply add the thermal con- 
tributions from the spin and charge sector. In that way, we can 
not discriminate the spin and charge velocities separately but 
their inverse sum, l/{l/vs + l/vc), also given in Table I. In 
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FIG. 4. (Color online) S2 as a function of block size for tlie half-filled 
chain at U /t — 6 and T/t = 0.06 using a grid of 51 = 4. The en- 
tropy data has been shifted such that S{Sl) = 0. The orange dashed 
line is a guide for the eye with a slope nTlA/{4:Vs) (see Equation 
(6)) that would correspond to purely thermodynamic entropy. Inset: 
S2/T for a half-system bipartition {I a ~ L/2) for different temper- 
atures from T/t = 0.04 to 0.08. The orange dashed line is a linear 
fit to the data at T/t = 0.06 from which we obtain the spin velocity 
to be Vs = 0.92 ± 0.05 (see also Table I). 



fact, we find good agreement between analytic results and our 
QMC calculations justifying our fitting ansatz. 



C. Mutual Information 

We calculate the mutual information I2 between two 
equally large blocks pinned at the left and right boundary. The 
blocks are grown simultaneously and have a size L/2 — A/2 
where A is the separation of the blocks (see Figure 6 for an 
illustration of the geometry). I2 for the quarter- filled chain 
at U/t = 4 is shown in Figure 6 for different system sizes 
and temperatures. First of all, I2 is sensitive towards the 2/cf 
oscillations that disappear for larger temperatures as it was ob- 
served already in Figure 3 for the Renyi entropies. As a func- 
tion of block separation, the mutual information is suppressed 
exponentially at finite T due to the presence of a finite correla- 
tion length and thus the correlations between two subsystems 
at, say T/t = 0.5, only extend over a few lattice sites. 

From Figure 5 one can also see that the mutual information 
at A = 0, henceforth denoted as e, increases with tempera- 
ture. A more detailed analysis, shown for the Mott insulator 
in Figure 6, reveals an interesting structure. Starting from the 
well known T — limit, e increases up to temperatures on 
the order of the superexchange scale, T k, J, where it ex- 
hibits a plateau. For T > 2 J, e shows a rapid increase again 
until it saturates eventually for temperatures much larger than 
U. This behavior seems counterintuitive at first sight because 
we expect that thermal fluctuations destroy correlations be- 
tween the subblocks since the correlation functions become 
more and more short ranged. In particular, the infinitely hot 
state is maximally mixed and the (reduced) density matrix 



is proportional to the unit matrix, p°°{A) = Z~^lz where 
Z — 4'-^* is the partition function and l^is d x d unit matrix. 
The systems considered here, however, are constrained by the 
total number of particles and zero net spin thus the particle 
number fluctuations of the two subblocks can not be totally 
uncorrected. The effect of particle number conservation on 
the complexity of infinitely hot states has been discussed in 
the context of matrix product operators^"-^' where it is found 
that representing p°° requires a finite matrix rank and the op- 
erator space entanglement entropy^- S^ of p°° (A) diverges 
logarithmically, i.e. S^ = 1/2 log i + const. In fact, we see 
that e for the unconstrained system - in Figure 6 we show data 
from grand-canonical simulations where the average particle 
number is one - is monotonically decreasing and tends to zero 
very quickly. 

The reduced density matrix for a subblock A of the T = 00 
state, p"^, and henceforth the mutual information of the in- 
finitely hot state denoted by e°° in the presence of the particle 
and spin constraints can be obtained by simple combinatorial 
arguments. For ease of simplicity, we focus on the case of a 
half system bipartition and unit filling but the results can be 
readily generalized to arbitrary fillings or partitions. First, we 
realize that the density matrix is a direct product of the spin 
components thus it suffices to consider one species of filling 
L/2 and we restrict p'^ to one spin sector In the second step, 
we write the p'^ as a sum over the number of particles con- 
tained in block A, p'^ — Z~^ ®„=o Pa^ where each sub- 
block is given as 



J") 



Pa =PL/2(i/2-n)V/2(")- 



(8) 



Vm (n) — (™) is the number of possibilities to put n particles 
onto m lattice sites and Z = p^ {L/2) is the partition function. 
The nth moment can then be calculated as 



tr[p; 
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[pL/2iL/2-n)rpL/2{n) 



(9) 



The Renyi entropy and thus e°° can then be evaluated using 
the definitions in Equations 2 and 7. The result is shown in 
the inset of Figure 6 comparing it to the T = result for 
U/t — 8. In particular, it is straightforward to show that e°° 
diverges like (a — l)/2 log L for large L where a is the Renyi 
index. 

Turning back to the results in the main panel of Figure 6 
where we normalized e with respect to e°°, we see that the 
plateau at T fa J corresponds to e w e°° /2. This value is 
expected if the spin degrees of freedom are completely disor- 
dered but the charge channel is still gapped out. Thus mea- 
suring e/e°° in the canonical ensemble can provide us with 
direct information about the relative energy scales in the prob- 
lem and whether we have some correlation apart from the 
constraints between the spin degrees of freedom. For large 
temperatures, U sets the energy scale for the saturation to 
e/e°° = 1. 

To probe this effects with cold atom experiments using the 
scheme based on a atomic gas microscope '^'^° one should note 
that the scheme proposed there does not allow a post selection 
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FIG. 5. (Color online) Mutual information I2 for a quarter-filled 
Hubbard chain at U/t = 4 and finite temperatures. The inset illus- 
trates the geometry of the two blocks that grow simultaneously with 
A sites between their edges, i.e. the block size is L/2 — A/2. The 
orange shaded area denotes the mutual information for a bipartition, 
e, that is detailed for the Mott insulator in Figure 6. 
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FIG. 6. (Color online) Mutual information for a bipartite splitting 
(A = 0) as a function of temperature for a Mott insulator at U/t = 8. 
The filled symbols correspond to data grand-canonical data with on 
average one particle per site. The data is normalized by Eoo (see text) 
and the vertical dashed line marks the superexchange scale J. Inset: 
e°° for a half-filled system as a function of system size. The black 
dashed line denotes the asymptotic scaling e°° — 1/2 log L + const, 
for large L. The green data points correspond to 25*2 for the ground 
state at U/t = 8. 



above, one has to ensure that the particle number fluctuations 
between pairs of tubes are small, e.g. as a consequence of 
the preparation process. For example one could prepare the 
state adiabatically, starting in a state with a pre-selected par- 
ticle number. The latter can be achieved e.g. with techniques 
realized in Ref. 73. This effect is also important in the context 
of particle number fluctuations, used for instance in measures 
like bipartite fluctuations'*^-^"*, where the constraint will also 
affect the finite-temperature behavior 



IV. DISCUSSION 

We used recent QMC methods to access the finite- 
temperature crossover of the one-dimensional Hubbard 
model. We can explore the domain up to very high temper- 
atures where the entropy becomes extensive. In particular, we 
accessed the second Renyi entropy to characterize quantum 
entanglement and thermodynamic dominated regimes by ex- 
plicitly comparing to ground state DMRG data. Characteristic 
features such as the 2fcF oscillations are suppressed rapidly 
and the purity is exponentially small both in temperature as 
in system size. Using the connection to the low-energy ther- 
modynamics of the Luttinger liquid description, the spin and 
charge velocities can be extracted from the thermal contribu- 
tion to the Renyi entropy. In view of future experiments, it 
is thus possible to see features of entanglement in the Renyi 
profiles for not too large systems and extract vital informa- 
tion about the underlying Luttinger liquid using a measure- 
ment scheme that is solely based on single-site resolved den- 
sity measurements"*. An important issue, that is beyond the 
scope of this investigation, is an analytical description of the 
thermal crossover regime of the Renyi entropies for finite sys- 
tems and we leave it for fuither studies. 

We also showed the mutual information for different block 
sizes that is also very short-ranged for intermediate tempera- 
tures and shows similar features as the plain Renyi entropies. 
In the presence of constraints, we also demonstrate that the bi- 
paitite mutual information can actually increase as a function 
of temperature and caiTies vital information about the nature 
of the fluctuations of the constrained and unconstrained sys- 
tem. 
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